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Abstract 



We propose a simultaneous solution to the strong CP problem and the 

SUSY phase problem based on parity symmetry realized when the supersym- 

_^ , metric standard model is embedded into a left -right symmetric framework at 

D \ a scale near 2 x 10^^ GeV, as suggested by neutrino masses and gauge coupling 



unification. In this class of models, owing to parity, SUSY contributions to the 
muon anomalous magnetic moment can be naturally large without conflict- 
^ ■ ing with the EDM of the electron and the neutron. The strong CP violation 

parameter 9 is zero at the tree level, also due to parity (P), but is induced 
due to P-violating effects below the unification scale. We estimate the in- 
duced ^ to be < 10^^^, if we adopt a constrained supersymmetric spectrum 
with universal scalar masses. In the more general SUSY breaking scenario, 
after imposing flavor changing constraints, we find ~ (10~^ — 10"^'^), which 
is compatible with, but not much below the present limit on neutron EDM. 
We also argue that potential non-perturbative corrections to 6 from quantum 
gravitational effects are not excessive in these models. 



I. INTRODUCTION 

One of the major problems of the Standard Model is a lack of understanding of the CP 
violating parameter 6 characterizing the QCD sector of the Lagrangian |]l[ . This parameter 
originates from the periodic vacuum structure of QCD and leads to an electric dipole moment 
(EDM) of the neutron, (i„ ~ (5 x 10^^^) 6 e-cm. The current experimental limit, dn < 
6.3 X 10~^^ e-cm |Q, then implies that 6 < 10"^^°. Why a fundamental parameter of the 
theory, 6, is so small compared to its natural value of order unity is the strong CP problem. 
A resolution of this problem is expected to provide important clues to the nature of new 
physics beyond the Standard Model. 

When the Standard Model is embedded into its minimal supersymmetric extension 
(MSSM) in order to solve the quadratic divergence problem associated with the Higgs boson 
mass, one runs into another CP problem - the SUSY phase problem, or the SUSY CP prob- 
lem. This problem owes its origin to the complex phases associated with the parameters 
in the soft SUSY-breaking sector of the theory - the /i term, gaugino masses, trilinear A 
terms, etc. Exchange of supersymmetric particles in loops will induce electric dipole mo- 
ments (EDMs) for the neutron (dn) and the electron (de) proportional to these SUSY phases. 
dn induced through gluino and squark exchange can be estimated to be (in ~ (1 x 10~^^) 
e-cm (300 GeV/MsusY)^sin0susY, where 0susy is a typical SUSY CP phase parameter and 
^susY is the gluino/squark mass. (For this estimate we used tan/3 = 5 and set the /i pa- 
rameter and the gluino mass equal to the squark mass.) In order to be compatible with the 
experimental limit on dn, the SUSY phase must obey 0susy < (5 x 10~^ — 5 x 10~^) for 
^susY = 300 GeV — 1 TeV, unless the gluino contribution is precisely canceled by some 
other diagrams [Q]. Why 0susy is so small, while the corresponding CP violating phase in 
the CKM matrix is of order one (in order to explain the observed CP violation in the K 
meson system), is the SUSY phase problem. 

The SUSY phase problem becomes more acute if we attribute the recently reported 2.6 
sigma discrepancy in the muon anomalous magnetic moment (a^) measurement Q] to su- 
persymmetry. Exchange of charginos/neutralinos and sleptons can account for the observed 
discrepancy in a^, provided that the masses of these particles are not more than about 
200 — 400 GeV. Now, if we replace the external muons in the diagram responsible for a^ by 
electrons, a large EDM for the electron will result, if the relevant SUSY phases are of order 
one. For example, if the SUSY contribution to a^ is 40 x 10~^°, de can be estimated by 
scaling of the lepton mass to be dg — (1.8 x 10~^^) sin^susY e-cm. The current experimental 
limit on de, viz., de < 4.3 x 10~^^ e-cm would require 0susy < 2 x 10~^. This bound 
will be tightened even further with the expected improvement in the limit on de by about a 
factor of 2 §. 

Solving the strong CP problem as well as the SUSY phase problem are therefore major 
challenges facing the (supersymmetric) Standard Model. 

Simple solutions to the strong CP problem can be found by postulating new symmetries 
of Nature, the most popular one being the Peccei-Quinn U{1)a symmetry 0]. Consistent 
implementation of this symmetry requires the existence of a new ultralight particle, the 
axion, which has eluded experimental searches so far. Combined laboratory, astrophysical 
and cosmological limits constrain the scale of Peccei-Quinn symmetry breaking to be in a 
narrow window, fa ~ (10^^ — 10^^) GeV [jl|. On the theoretical side, global symmetries such 



as U{1)a have come under suspicion since it is believed that nonperturbative quantum grav- 
itational effects will violate all global symmetries. If so, quantum gravity would destabilize 
the axion solution [|^ unless one allows for extreme fine-tuning of parameters - the very 
problem one set out to avoid in postulating the new symmetry.^ In any case, this solution 
has nothing to offer to the second CP problem, the SUSY phase problem. 

An alternative to the axion solution to the strong CP problem is parity invariance 
realized at a momentum scale vr much above the weak scale.0 If the Standard Model is 
embedded into a left-right symmetric gauge structure at vr, parity (P) invariance can be 
consistently imposed on the Lagrangian. In this case, the 6 term in the QCD Lagrangian, 
g'^/{32'ir'^)6GG, will be zero since it violates parity. The physical parameter 6 involves also 
the phase of the fermionic determinant, and is given by 

9 = e + ArgjDet (M^Mrf)} - 3Arg(M3) . (1) 

Here Mu^d are the up-quark and the down-quark mass matrices and Mg is the gluino mass. 
Owing to parity invariance, the matrices M„ and M^ become hermitian, and the gluino mass 
becomes real. As a result, ^ = at tree level in this class of models. These models would 
thus have the potential to solve the strong CP problem. The fact that low energy weak 
interactions do not respect parity symmetry means that one must do additional work to see 
if 6 induced through quantum effects is sufficiently small. As we shall explicitly demonstrate 
in this paper, this is often the case. 

The purpose of this note is to provide a simultaneous solution to the strong CP problem 
and the SUSY phase problem using parity symmetry. We shall demonstrate by explicit 
calculation that the induced 6 in these models is well within the experimental limit, if a 
constrained supersymmetric spectrum is adopted with universal squark masses and propor- 
tional A terms. Even in the more general scenario for supersymmetry breaking, we shall 
see that 6 is in the acceptable range of ~ 10"^'' — 10~*, after imposing flavor changing con- 
straints. Parity invariance also makes the phases of the SUSY breaking parameters naturally 
small. These phase parameters are zero at the scale vr and their induced values at the weak 
scale through quantum corrections are well within the experimental limits arising from the 
neutron and the electron EDM. Thus this class of models can naturally explain the observed 
discrepancy in the muon anomalous magnetic moment a^, without inducing unacceptably 
large EDM for the neutron and the electron. 

In the class of models presented here the scale of parity restoration is in the range 
vr ~ 10^^ — 10^^ GeV. This is close to the grand unification scale where the three gauge 
couplings of the Standard Model are observed to unify in a supersymmetric context. The 



^For example, quantum gravity can induce a dimension 5 operator in the scalar potential 5|<1>|^($+ 
$*)/Mpi, where g is a dimensionless coupling, and ^ is the scalar field responsible for U{1)a 
symmetry breaking. Such a term will induce a nonzero 9 given by ^ ~ gf^/{AQ(-,j-)Mpi) ~ 10 g. 
The resulting constraint on the coupling g from neutron EDM is quite severe: g < 10"^'^. 



^For models with CP invariance, see [10|. Related applications in SUSY context has been discussed 
inRef. Q. 



left-right symmetric gauge structure and the numerical value of the scale vr are indepen- 
dently suggested by experimental evidence for neutrino masses: Small neutrino masses aris- 



ing through the seesaw mechanism |]T2| is natural in the left-right framework, and the scale 
vr is consistent with the inferred value of Uj. mass from atmospheric neutrino oscillations. 
Parity symmetry as a possible solution to the strong CP problem in supersymmetric 



contexts has been studied in earlier papers ||T3[ . The models presented here are significantly 



improved versions in this regard. In earlier works |TB|, it was found that a consistent solution 
to the strong CP problem required the scale of parity restoration to be in the multi-TeV 
range, which would appear to not go well with neutrino masses and gauge coupling uni- 
fication, unlike in the models presented here. As we shall see, potential non-perturbative 
corrections to 6 from quantum gravity are under control here, even with vn near the unifica- 
tion scale. Parity as a solution to the SUSY phase problem has also been studied in earlier 
papers [|14| , |T5| , where it has been shown that the EDMs of the neutron and the electron 
remain very small in the parameter space that fits e and e' in the Kaon system. Unlike in 
these earlier works which preferred a non-KM mechanism for the Kaon CP violation, the 
models here allow for the conventional KM CP violation. This is facilitated by a novel real- 
ization of the doublet-doublet splitting - the mechanism that makes one pair of Higgs(ino) 
doublets light and all other pairs heavy, so that at low scale the spectrum of the theory is 
identical to that of the MSSM. Major differences of our models compared to the MSSM are 
that here (i) SUSY phases are naturally small, (ii) the strong CP problem is absent, and 
(iii) small neutrino masses are naturally present. 

II. BASIC OUTLINE OF THE MODEL 

The basic framework of our model involves the embedding of the MSSM into a minimal 
SUSY left-right gauge structure at a scale vr close to the GUT scale. The electroweak gauge 
group of the model is SU{2)l x S'f/(2)/j x U{1)b^l with the standard assignment of quarks 
and leptons - left-handed quarks and leptons {Q, L) transform as doublets of SU{2)i, while 
the conjugate right-handed ones {Q'^, L'^) are doublets of SU{2)fi. The quarks Q transform 
under the gauge group as (2, 1, 1/3) and Q'^ as (1,2, —1/3), while the lepton fields L and 
L'^ transform as (2,1,-1) and (1,2,-M) respectively. The Dirac masses of fermions arise 
through their Yukawa couplings to one or more Higgs bidoublet $(2,2,0). The SU{2)r x 
U{1)b-l symmetry is broken down to f/(l)y in the supersymmetric limit hj B — L = ±1 
doublet scalar fields, the right-handed doublet denoted by ^'^(l^S, —1) accompanied by its 
left-handed partner x(2, 1, 1). Anomaly cancellation requires the presence of their charge 
conjugate fields as well, denoted as x'^(l;2, 1) and x(2, 1,— 1). The vacuum expectation 
values (VEVs) {x'^) = (x^) = vr break the left-right symmetry group down to the MSSM 
gauge symmetry. 

This embedding of MSSM into a left-right framework provides a simple solution to the 
SUSY phase problem. To see this, let us note the transformation of various fields under 
parity symmetry: Q ^ Q"*, L ^ L"*, $ ^ $^, x ^ X"* ■, X ^ X"*-, G ^ G*, B ^ B*, 
Wl •«-»• W^ and 9 ^^ 9. Here (G, B, Wl,r) are the vector superfields associated with SU{3)c, 
B — L and SU{2)lji respectively, 9 is the fermionic variable, and the transformation applies 
to the matter superfields as a whole. Invariance under P makes the Dirac Yukawa couplings 
of quarks and leptons and the associated SUSY breaking trilinear A-terms hermitian. The 



gluino and the B — L gaugino masses become real, the mass of the bidoublet field $ as well 
as the corresponding bilinear S/i term are real, and M^ = M^, Mi^^ being the masses of 
the SU{2)l^ji gauginos. This resolves the SUSY phase problem, since all the relevant SUSY 
phases are zero at the scale of parity restoration [Q . Renormalization group extrapolation 



induces very small phases in the SUSY breaking parameters of the MSSM, but these in- 



duced values are well consistent with experimental limits arising from dn and d^ ||T5[. It is 
worth noting that this solution to the SUSY phase problem will be valid even in a general 
context of SUSY breaking, for example, without assuming universality of scalar masses and 
proportionality between the A terms and the Yukawa couplings. Potential contributions to 
the EDM of the neutron and the electron will be proportional to the diagonal entries of 
the respective A matrix or the squark/slepton mass-squared matrix. Both matrices being 
hermitian, these contributions vanish above the scale vr. Since the gaugino masses are all 
real (assuming gaugino mass unification that occurs in various scenarios of SUSY breaking 
even without a unifying group, or if the left-right gauge theory is embedded into a higher 
symmetry group such as [SUi?))]^, 5*0(10) or Eq), dn and d^ proportional to the phases of 
the gaugino masses also vanish above vr. 

A. The doublet doublet splitting mechanism 

In order to make the supersymmetric left-right gauge theory fully realistic, a mechanism 
should be found that keeps one pair of Higgs doublets light at the weak scale (to be identified 
with the Hu and Hd fields of MSSM) and any remaining pairs of Higgs doublets superheavy 
at the scale vr. The simplest possibility would appear to be to introduce just a single Higgs 
bidoublet $(2, 2, 0) which gets a mass of order SUSY breaking scale. However, this is not 
the minimal scenario from the effective low energy point of view, since in that case the up 
and down quark mass matrices become proportional, leading to vanishing quark mixings at 
tree level.0 

In this subsection we shall present a new mechanism for doublet-doublet splitting. It 
involves two bidoublet Higgs fields and is achieved without any fine-tuning of parameters. 
The solution to the SUSY phase problem is preserved based on parity symmetry alone. 



^Realistic quark mixings can be induced through the gluino and chargino loops, provided that 



the trihnear A terms have a flavor structure different from that of the Yukawa couplings |15|,16|. 



Consistency with flavor changing processes would require that tan/3 be not too large |15|, tan/3 < 
10, which excludes the simplest scenario where tan/? = rat/rrn, ~ 60. Values of tan/? smaller than 
rrit/rrn, may be obtained if the $ fleld mixes with some other superheavy doublets of the theory 
in a parity violating manner (if such mixings conserve P, tan/? = mt/mb will prevail). In this 
case the effective fi and B/j, terms are potentially complex, which would spoil the solution to the 
SUSY phase problem based on parity. To maintain this solution, in earlier works |0] we assumed 
invariance under charge conjugation symmetry C, in addition to P, which allows for reality of the 
effective fj, and Bfj, terms. In such a scenario the CKM phase will be zero due to C invariance and 
the observed CP violation in the Kaon system is explained through supersymmetric gluino/squark 
diagrams. 



As we shall elaborate further in the next subsection, in this scenario the quark mixings 
arise naturally at tree level, tan/3 can be smaller than mt/mb, symmetry breaking occurs 
at the renormalizable level without any pseudo-Goldstone bosons, and neutrino masses are 
correctly reproduced with v^ near the unification scale. Furthermore, the effective low 
energy theory is just the MSSM, with a natural understanding of the weak scale value of 
the fi parameter which remains real. 

Consider the following form of the superpotential involving two bidoublet field ^a, ci = 
1, 2 and the left-handed (x + x) ^^^ and the right-handed {x'^ + x'^) doublets of the theory: 

W = XaX^aX' + Kx^aX' + Mxx + M^x't ■ (2) 

Parity invariance makes the couplings Aq and A'^ real, since $ -^ <l>^ and x "^ X^* under P. 
The mass term M and the VEV (x'^) = (x^) = ^R ^i^re complex in general. After SU{2)ji 
breaking, this superpotential leads to a mass matrix for the doublets given in Eq. (3). We use 
a notation in which the rows denote (0mi, (pu2, x) fields and the columns denote {4>di, (pd2, x) 
fields where 0„i and 0^1 are the up and down type Higgs doublets from $i, etc. 



M, 



DD 
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\2Vr I . (3) 

\X'iVr X'^Vr 




This mass matrix leaves one pair of Higgs doublets massless, while giving mass of order 
vr to the second pair. Since Aa and A^ are real, the effective /x term of the hght doublets 
becomes real. To see this, observe that the low energy MSSM doublets are given by H^ = 
cos a„ (pui + sin a„ 0^2 and Hd = cos aa (pdi + sin ad 4>d2, where tan a„ = A1/A2 and tan ad = 
A'^/Ag. Note that H^ and Hd are real linear combinations of $„, which helps in inducing a 
real /x. The superpotential of Eq. (2) by itself does not lead to a /i term, which gets induced 
only after SUSY breaking. There are two sources that induce the /z term: 

1. Kahler potential terms of the form X'^f^Jd'^O^—Tr^^a^b), where Z is a gauge singlet 
whose Fz ^ breaks supersymmetry. We assume that Z is parity even, i.e., Z —>■ Z* 
under P. The coupling matrix A" is therefore hermitian. After supersymmetry breaking 
this term will lead to a real /x-term, as desired. This also provides a reason why the 
/i-term is of order of the electroweak scale |T^. The Bfi term arises from the term 
J d'^O^TT Tr($a$b), which is also real due to parity. 

2. A second mechanism that generates real /i and Bfi terms makes use of the superpoten- 
tial couplings involving a visible sector singlet S": W D K,S{e^^x'^X'^+^~^^XX~M'^)- Such 
a coupling can break the left-right gauge symmetry down to the MSSM symmetry at 
the renormalizable level without leaving any pseudo-Goldstone bosons. Owing to par- 
ity, under which S —* S*, the parameters k, M^ are real in this superpotential coupling. 
The field S also has the following coupling to the bidoublets $„: W D /iabTr($a$fc)S'. 
In the SUSY limit, (S) = 0, (x^) = (x'^) = "^-R; {x) = {x) = 0, which breaks parity 
spontaneously. S pairs up with the neutral component (x'^ + x'^)/"\/2 to form a multi- 
plet that has mass \P1kvr. After SUSY breaking, the coupling J d'^9^^^ — will induce 
a tadpole in Re(S') scalar of order f|,m3/2- A VEV (Re(S')) ~ ^-3/2 will result, that 
provides a real fi term for the bidoublet fields. It is crucial to note that by redefining 



the y^ filed, the coefficient of the tadpole in S can be made real without introducing 
any phases elsewhere (see Eq. (4) below). If the imaginary component of S had also 
a tadpole, (Im(S')) 7^ 0, which will lead to an effective complex /x term. 

B. The full Lagrangian 

Now let us implement the doublet-doublet splitting mechanism just described. We shall 
see that there is a discrete anomaly-free Z4 R-symmetry that achieves this goal within a 
minimal version of the left-right model (viz., using two bidoublets <l>a, one left-handed (x) 
and one right-handed (x'^) SU{2) doublets along with their conjugate (x+X'^) ^^^ the singlet 
S). Their transformations under P has been given earlier, with S -^ S*. All the desired 
terms, including the Majorana mass terms for the right-handed neutrinos are allowed by 
this Z4, and the unwanted terms that can potentially make the magnitude of fi too large, or 
induce excessive CP phases to upset the strong CP and the SUSY phase solutions will be 
prevented. The Z4 is broken at the scale vr, but a Z2 remnant remains, which is identified 
as the usual R parity of the MSSM. This Z2 will guarantee the stability of the proton. 

Under the Z4 R symmetry, the superpotential changes sign {W -^ —W), as do d'^6 and 
d^6. The gaugino fields transform as Xa -^ — Aa, quarks and leptons are even, <l>a : —1, x '■ 

h X ■ -h X"" ■ -h X'' ■ +h S : -1. 

The gauge invariant superpotential consistent with this Z4 R symmetry is 

W = KQ^aQ' + KL^aL'' + Kx^aX' + Kx^aX" + U LLxx + f* L' L' x' x') / Mpi + 

f^Sie^'x" + e-'^XX + aS'' - M^) + /i,,Tr($,$5)5 . (4) 

This superpotential induces tree level CKM mixings since the light MSSM doublets Hu,d 
are parity asymmetric linear combinations of the two bidoublets. The / couplings give rise 
to Majorana masses for v^ of order v^/M-pi. For vr ~ 10^'* — 10^^ GeV, the magnitude of 
the light neutrino masses are in the right range to explain the atmospheric and the solar 
neutrino oscillation data. / could have its origin in quantum gravity, but it could also arise 
from integrating out singlets which have masses of order Mpi, e.g., through {LN x + L'^ N^x'^) 
couplings where {N,N^) are the singlets with Z4 charges {i,—i). Their Majorana masses 
[A^^ + A^'^ ] preserve the Z4 symmetry. 

In the SUSY limit, we have (S) = 0, (x'^) = (x*^) = ^ with all other fields having zero 
VEVs. As noted earlier, after SUSY breaking, the real component of S gets an induced 
VEV of order 1713/2. (Note that the phase in the x'^X'^S scalar coupling can be made real 
by redefining x'^ field. This redefinition does not induce new phases anywhere else. The 
tadpole in S that is induced after SUSY breaking is therefore real, making only Re(S') to 
be nonzero.) That gives a hermitian fiab terms through the last couplings of Eq. (4), or to 
real /i parameter. We can also have the coupling X'^f^J d'^9 ^ '^^^ '''' where Z is the spurion 
field that breaks SUSY. This also leads to real fi term of the right order of magnitude. 
(Note that Z is parity even, Z -^ Z* under P. Fz is then expected to be real, which would 
leave parity unbroken. For example, in the Polonyi model of hidden sector SUSY breaking, 
W D fi^{Z + P), where fi^ is real due to parity. Fz = jj? is therefore real. We anticipate the 
reality to Fz to hold even in a more general scenario for SUSY breaking.) 

The superpotential of Eq. (4) reproduces the doublet-doublet mixing matrix of Eq. (3). 
In Eq. (3), the (3,3) entry is of order 7723/2 now, being proportional to {S). It does not 

7 



correspond to any new particle having mass of order ^3/2, since x pairs with the heavy 
doublet in $„ and has a mass of order vr. 

We shall now show that the Z4 R symmetry is an anomaly free discrete gauge symmetry 



18| . This makes it aesthetically more pleasing, as it may have its origin in a true gauge 
symmetry. It also protects the Lagrangian from receiving uncontrollable quantum gravita- 
tional correction. To see the anomaly freedom, let us assume that the Z4 arose from a true 
gauge U{1)r symmetry. The U{1)r then should be anomaly free. If the U{1)r symmetry is 
broken at the Planck scale by scalar VEVs that are integer multiple of ±4 (upto an overall 
normalization factor), a residual Z4 symmetry will survive. The U{1)r anomaly cancellation 
will in general require introduction of additional fermionic fields. The crucial question is 
then if these extra fields can be removed from the low energy spectrum by giving them Z4 
invariant masses. To address this, let us embed the Z4 into a U{1)r in the obvious fashion, 
by assigning U{1)r charges as follows: (We display the R charge of the superfield, which 
is the same for the scalar component, but the fermionic component will have its R charge 
shifted by —1.) Quarks and leptons: 0, x '■ +1; X'^ '■ ~^i X '■ ~1; X'^ • +1; ^a '■ +2, S : +2. 
The superpotential W has R charge of +2, and the gauginos have R charge of +1. 
With this assignment, one can compute all the mixed anomaly coefficients: 



SUi-ifc X U{1)r 
SU{2)1 X Uil)R 
SU{2)1 X U{1)r 
U{1)1_^ X U{1)r : -3 - 2Ng (5) 



?>-2Ng 
3-2Ng 

3-2Ng 



Here Ng = 3 is the number of generations. The 3 in first term arises from gluino loop. The 
3 in second term is from Wino (+2) and Higgsino (+2 from the two bidoublets and —1 from 



x). We have used the conventional 5*0(10) B — L normalization, J3/2(B — L)/2 being the 
normalized generator. 

Since all the nonabelian mixed anomaly coefficients are equal, we can try to cancel them 
by the Green-Schwarz mechanism [^]. We shall make the abelian mixed anomaly coefficient 
to be also equal to facilitate this. This can be achieved by adding a pair of singlets which 
are (2,3) + (—2,3) under {B — L,R). Then the U{l)%_i^ x U{1)r anomaly also becomes 
(3 — 2Ng). A mass term for these singlet fields will have R charge of +6 (scalar component), 
so it must be accompanied by a Higgs field with charge —4. That breaks U{1)r to Z4, as 
desired. 

Finally, there is t/(l)|.x U{1)b-l anomaly, which has a coefficient — 8-y/3/2 in this model. 



We use the overall R-normalization (level of U{1)r) to make this equal to the other anomaly 



coefficients. This normalization factor is then found to be (1/4) w 3/2 - very similar to B — L 
normalization. The [f/(l)/j]^ anomaly also has the same coefficient (equal to —3 for A^^ = 3) 
if some singlet fields contribute +24 in the cubic anomaly. 3 singlets with fermionic R- 
charge of +2 will do this job. Mass terms for these singlets will carry R-charge of +6 (scalar 
component), so to make it +2, we must multiply by a Higgs with R-charge of —4. Again, 
we see that the Z4 is left unbroken by this Higgs. This way of canceling anomalies is not 
unique, but to establish that Z4 is discrete anomaly free, any one example would suffice. 



III. CALCULATION OF INDUCED 9 

As noted earlier, parity invariance implies that the QCD Lagrangian parameter ^ = 0, the 
gluino mass is real and the quark mass matrices Mu^d are hermitian at tree level. Therefore 
^ = at tree level. Since parity is broken at f r, a nonzero value of 9 will be induced at the 
weak scale through renormalization group extrapolation below vr. We shall estimate this 
induced 9 in two scenarios for SUSY breaking. The first is the constrained MSSM scenario 
where the squark masses are degenerate at the unification scale with the trilinear A matrices 
and the corresponding Yukawa coupling matrices being proportional. The second scenario 
has a more general SUSY breaking spectrum without universality or proportionality, but 
the experimental constraints arising from flavor changing processes will be imposed. We 
flrst turn to the correction to 9 arising from the non-hermiticity of the Yukawa coupling 
matrices which applies to both these scenarios. 

1. 59 from non— hermiticity of the Yukawa coupling matrices: 

At the scale fR, the up and down Yukawa coupling matrices are hermitian owing to 
parity. The VEVs of the MSSM flelds Hud are real since the B^ term is real, so that the 
quark mass matrices M^ and Md are hermitian at vr. They will develop non-hermitian 
components at the weak scale, owing to renormalization group evolution below vr. The 
induced 9 will have the general structure given by 

59 = ImTr[AM„M~^ + AMaM^^] - 3lm{AM~gM^^) (6) 

where Mu^d,g denote the tree level contribution to the up-quark matrix, down-quark matrix 
and the gluino mass respectively, and AMu^d,g are the loop corrections. To estimate the 
corrections from AMu and AMd, we note that the beta function for the evolution of Y^ 
below Vr is given by (3yu = ^M/(16vr^)(3yJ"yu + Yj^d + Gu) with the corresponding one for 
Yd obtained by the interchange F„ ^^ Yd and Gu -^ Gd- Here Gu is a family-independent 
contribution arising from gauge bosons and the Tr{Yj^Yu) term. The SY^Yu term and the G„ 
term cannot induce non-hermiticity in Yu, given that Yu is hermitian at vr. The interplay 
of Yd with Yu will however induce deviations from hermiticity. Repeated iteration of the 
solution with Yu oc YuYjYd and Yd oc YdY^Yu in these equations will generate the following 
structure: 

where Mu is the unification scale. Here ci and C2 are order one coefficients which are not 
equal since the fiavor independent parts G„ and Gd are not the same for the evolution 
of Yu and Yd (hypercharge gauge couplings and the tau lepton Yukawa couplings differ- 
entiate the two.) These contributions to 59 are very high order in the Yukawa couplings 
since the trace of products of two hermitian matrices, having the form Im Tt{YuYPY^Y^..) 
contains an imaginary piece only at this order. To estimate the induced 9, we choose a 
basis where Yu is diagonal, Yu = D and Yd = VD'V^ where Z)„f„ = diag(?72„, rric, mt), 
DdVd = diag{md, rris, nih) with V being the CKM matrix. The Trace of the first term 
in Eq. (7) is then Ira^Df DfDj D[ VijVkiV*iVj^j) . The leading contribution in this sum is 



{rn^rnl'mfrri'l) / {vlv^)lin{VcbVtsV*gVfl) . The second Trace in Eq.(7) is identical, except that 
it has an opposite sign. Numerically then, 

66 ^3x 10"2^(tan/?)^(ci - Cs) (8) 

where we have used the running quark masses at rrit to be {mi;,mc,mb,ms) = 
(166, 0.6, 2.8, .063) GeV. Clearly, 69 is very small, even for tan/5 = 50 its value is 10~^^, 
much below the experimental limit of 10~^° from neutron EDM. 

Since the up(down)-quark mass matrix is a product of Yu(Yd) and the VEV Vu{vd), the 
mass matrix can become complex if the VEV Vu{vd) is complex. If the bilinear soft SUSY 
breaking parameter Bfi becomes complex in the process of evolution below vr, this will 
happen. By analyzing the RGE for the Bn parameter, one sees that it involves traces of 
(Y^Yy^) and (YjYd) or their products - in the case of universal squark masses and proportional 
A terms [A^ oc Y^, A^ oc Y^). We are again left with two hermitian matrices [Y^, Yd), with 
all other effective parameters being real. The imaginary component of the trace that induces 
a phase in B^ is then given at lowest order by an expression analogous to Eq. (7). The 
estimate on 66 is of the same order as before, 66 ~ 10~^^(tan/3)^. 

2. 66 from finite correction to the quark and the gluino masses: 

To compute 66 arising from the finite corrections to the quark mass matrices and the 
gluino mass (which are not contained in the RGE evolution), we must specify the SUSY 
breaking spectrum. The simplest approximation is to assume universality of scalar masses 
and proportionality of A-terms and the respective Yukawa couplings at the Planck scale. 
This can be justified in models such as the ones with gauge mediated supersymmetry break- 
ing [^. In this case, the whole theory at the weak scale is characterized by only two Yukawa 
coupling matrices Yu^d- Furthermore, all other MSSM parameters are real in the effective 
low energy theory below vr- Because of this property it is very easy to estimate the lowest 
order contribution to nonvanishing 6 in terms of the coupling matrices. 

Consider the finite one loop corrections to the quark mass matrices. A typical diagram 
involving the exchange of squarks and gluino is shown in Fig. 1. There are analogous 
chargino diagrams as well. In Fig. 1, the crosses on the Q and Q^ lines represent (LL) and 
(RR) mass insertions that will be induced in the process of RGE evolution. From this figure 
we can estimate the form for AMu = "^rrAAurn?^^ where Q is the squark doublet and u'^ is 
the right-handed singlet up squark. Without RGE effects, the trace of this term will be real, 
and will not contribute to 6. Looking at the RGE for m\ upto two loop order, we see that 
for the case of proportionality of A^ and Y^., m\ gets corrections having the form mlY^ or 
mlY^ or mlYuY^Yu- Therefore in AMuM~^, the M~^ always cancels and we are left with 
a product of matrices of the form Y^Y^Y^Y^ ■ ■■. A similar comment applies when we look 
at the RGE corrections for m^ or Au- If the product is hermitian, then its trace is real. So 
to get a nonvanishing contribution to theta, we have to find the lowest order product of Y^ 
and Y^ that is non-hermitianQ and we get 



■^Similar reasoning was used in the standard model and super symmetric models in earlier papers 
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2a, nn{Mu/M w] 



S0 = -^ (^ ":^i;"^ ) ( fcilmTr [r„^F,^F.^yi] + fc.ImTr [Y.^Y^Y.X^] ) (9) 



where /ci 2 are calculable constants. The numerical estimate of this contribution parallels 
that of the previous discussions, 66 ~ {ki — ^2) x 10~^^(tan/3)^. The contributions from the 
up-quark and down quark matrices tend to cancel, but since the d'^ and the u'^ squarks are 
not degenerate, ki 7^ /c2 and the cancellation is incomplete. 



O^- <V--^^Q^ 
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Q/ \Q^ 
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FIG. 1. One-loop gluino/squark exchange diagram contribution to the quark mass matrix. 
The crosses on the scalar lines correspond to mass insertions. 

In Fig. 2 we have displayed the one-loop contribution to the gluino mass arising from the 
quark mass matrix. Here again one encounters the imaginary trace of two hermitian matrices 
Yu and 1^, in the case of universality and proportionality of SUSY breaking parameters. Our 
estimate for bQ is similar to that of the quark mass matrix of Eq. (9). 
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FIG. 2. One loop diagram that inuduces a phase in the gluino mass. 

3. Induced Q with general SUSY breaking terms: 

In this subsection, we study the more general SUSY breaking scenario where soft SUSY 
breaking terms involving squarks is given by 

l^ssB = E ^^ml(P + QA^HuU^ + QAdHj^ + h.c. (10) 

For the model under study, at the vr scale, the constraint is that m^ = m^-^ = m|^ = m^^ 
due to parity invariance. A^^d are arbitrary hermitian matrices, and the squark mass matrices 
can have non-trivial flavor structure. 

11 



In this case, the lowest order correction to 59 from one loop contributions to quark masses 
(Fig. 1) is given by 

for / = u^d. This contribution vanishes since the matrix m^Afm^ and Yf^ are both 
hermitian. The next leading contribution has the form 

This contribution arises from Fig. 1 by inserting «^^y„^ arising from the RGE equations in 
one of the squark lines. Since this trace involves three arbitrary hermitian matrices, it is 
not real in general. To estimate this contribution, we have to make some assumption about 
the non-universality in rn^ and the non-proportionality in A and the Yukawa coupling 
matrix. As for the A term, the most natural choice will be to assume that it has the 
same hierarchical structure as the Yukawa couplings. Such a form would be suggested 
by flavor symmetries. Thus, we shall take A23 ~ £^33, where e is a small parameter, of 
order Vcb ~ 1/30. Such a choice will guarantee that there is no excessive FCNC processes 
mediated by squarks. As for the squark mass matrices m^, we take it to be approximately 
proportional to a unit matrix, with correction terms that are not large. This is as suggested 
by nonabelian horizontal symmetries |Q. The leading contribution from Eq. (12) to 6 
arises when we use index (3,2) for the first m|, (2,3) for A^ and (3,3) for the rest. Using 
v^k/mo ~ 1/5, As3/mo ~ 1/10, A23 ~ lO'^Aga, (m|)23 ~ lO'^m^, we find 59 ~ lO'^. This 
is a conservative estimate and yet it is encouraging that we are close to the present upper 
limit on 59 of 10~^ to 10~^°. To be completely consistent with the neutron EDM limit, 
we should have the relevant phase to be of order 0.1, or the off-diagonal entries somewhat 
smaller than allowed by FCNC constraints. Since such departures from natural values need 
be only mild, we feel that this scenario is also quite viable. It is interesting that in this 
scheme, dn is not much below the present experimental limit, while dg is well below the 
current limit. 

Let us now address the contribution to 9 arising from the induced phase of 
the gluino mass. The leading contribution (see Fig. 2) in this case is given 

— 2 

by 59 :i^ ^^lmTT{AuYu) ^f^. ■ Without RGE running, this trace is real since 
Yu and A^ are hermitian. Allowing for RGE running, we estimate 59 ~ 
^ [ln(Mc//M^)/(167r2)]ImTr(M?y,A„)^. Taking the (2,3) entries of M? and A„ to 
be 10~^ times that of the respective (3,3) entries, and with A33 = ttiq/IO, we arrive at 
59 ~ 10~^ — 10~^° for Vwk/iTT'O ~ 1/5. This is again not far from the present upper limit and 
with a mild fine-tuning of parameters, of order 10%, one gets the desired solution to the 
strong CP problem. 



IV. PLANCK SCALE CORRECTIONS 

One interesting aspect of the model presented here is that it is quite safe from poten- 
tially large corrections to 9 induced by quantum gravity. If it is assumed that the high 
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scale parity conserving theory originates from a more fundamental theory, one can expect 
nonrenormalizable operators in the theory suppressed by the mass scale associated with the 
fundamental theory. Such corrections to 6 will respect the gauge symmetry as well as the 
anomaly free Z^ discrete gauge symmetry. We should ensure two things: (i) The effective /i 
term induced by quantum gravity is not more than the weak scale, and (ii) The quantum 
gravity induced phases which may not respect parity do not upset the solution to the strong 
CP problem. All other constraints, such as the solution of the SUSY phase problem, will be 
automatically satisfied once these two are taken care of. 

As for the magnitude and the phase of the effective fi term, the most relevant higher di- 
mensional operator suppressed by Planck mass that is invariant under the gauge symmetries 
and the Z^ symmetry isW D Kab'TT^{^a^b)x'^X'^S/Mpi. The magnitude of the resulting /i term 
is KabVJiMsvsY / Mpi ~ 10~^MsusY- Clearly, this is very small correction to the magnitude 
of /i. Suppose that quantum gravity does not respect parity symmetry. The coefficients Hab 
will then be non-hermitian. The phase of the /i term will then be arg(/i) ~ 10~*. Through 
the gluino diagram this will lead to ^ ~ 10"^'', which is consistent with dn limit. This shows 
that the complex couplings Kab can be of order one. 

The quark mass matrices can also have corrections from Planck scale physics. The most 
relevant term is W D QQ'^^x'^X'^ / ^pi^ which will induce corrections of order 10~^f^fc for 
some of the quark masses. The matrix structure need not be hermitian if quantum gravity 
violates parity. We suspect that gauged flavor symmetries (discrete or continuous) must 
exist in the underlying theory, or else the light fermion masses can become too large from 
quantum gravity. Most likely the estimate of lO"^^^;^ will apply for the third generation. If 
the coefficient of this non-renomralizable operator is of order 10~^ — 10~^, the solution to 
the strong CP problem via parity will be preserved. Note that the superpotential coupling 
W D QxQ'^X'^/Mpi is not invariant under the discrete Z^, unless accompanied by another 
factor S/Mp\. The correction to the quark mass matrix from this term is extremely small 
AM„ ~ {Msusv/Mpifv^k- 

We have verified that all other Planck induced corrections are much below the experi- 
mental limits on 6. 

Question can be raised as to the form of SUSY breaking parameters and if they indeed will 
respect parity symmetry. A complete answer to this will have to await a full understanding of 
non-perturbative SUSY breaking which is lacking at the moment. We note that perturbative 
gravity, which is utilized in conventional supergravity models of SUSY breaking may well 
respect parity - we have given an example in Polonyi model. A second example is gauge 
mediated SUSY breaking. If SUSY is broken at a scale of 10"^ — 10^ GeV, quantum gravity 
corrections for the fi term and the A term, which will be of order gravitino mass, will of order 
10^^° — 10"*^ GeV. Even if they are complex and non-hermitian, the strong CP problem will 
be solved, as the induced 6 will be of order 10~^° — 10~^^. We may use one of the other 



proposed solutions to generate a /i term of the weak scale in this case |^ . If the messenger 
fields do not couple to the fields x'^) X'^y they will not feel the effects of parity breaking, 
although parity is broken at f^ ~ 10^^ GeV. The effective SUSY breaking parameters will 
then obey the constraints of parity. 

The d = 5 baryon number violating operator QQQL/Mpi in the superpotential is for- 
bidden in this model by Z4 symmetry, but the operator QQQLS/Mp^ is allowed. If the 
associated couplings are order one for the light generations, we estimate proton lifetime 
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induced by these operators to be Tp ~ 10^'' yr. 



V. CONCLUSIONS 

We have shown in this paper that it is possible to embed the supersymmetric Standard 
Model into a parity-symmetric framework at a unification scale of 2 x 10^^ GeV in a simple 
way. Such an extension is well motivated by the data on neutrino oscillations as well as gauge 
coupling unification. We have demonstrated that this embedding can naturally solve the 
strong CP problem and the SUSY phase problem simultaneously. The effective low energy 
theory is the MSSM, but with naturally small phases for the SUSY breaking parameters 
along with order one phase in the CKM matrix. Thus it allows for large SUSY contribu- 
tions to the muon g — 2, as indicated by experiment, without violating the bounds on the 
electron EDM. The induced 6 in these models depends strongly on the way SUSY breaking 
is communicated. With universality of squark masses and proportionality of the A terms, 
we found 6 < 10~^^, while with maximal deviation from universality and proportionality 
consistent with FCNC constraints 6 ~ 10^^° — 10^^. In the latter case, the neutron EDM 
should be soon accessible, while de will be much smaller than the present experimental limit. 
We have also shown that potential corrections induced by quantum gravity are under con- 
trol in this class of models. Since left-right gauge symmetry is realized at a scale v^ ~ 10^^ 
GeV, evolution of couplings between vr and Mpi can induce flavor changing neutral current 
processes which are in the interesting range for current and future experiments. We plan to 
study this issue in detail in a forthcoming publication. 
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